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and the so-called 1 + 3 gauge-invariant approach to cosmological perturba- 
tions. We restrict our considerations to linear perturbations, but our intent 
is to set the stage for the extension to second order perturbations. 
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1 Introduction 

Currently, increasingly accurate observations are driving theoretical cosmology to- 
wards more sophisticated models of matter and the study of possible nonlinear 
deviations from FL cosmologyQ Motivated by this state of affairs, in a recent paper 
(Uggla and Wainwright (2011), hereafter referred to as UW), we initiated a program 
of research whose long term goal is to provide a general but concise description of 
nonlinear perturbations of FL cosmologies that will reveal the structure of the gov- 
erning equations, and hence facilitate their analysis. In furthering this goal one is 
faced with making three choices. First, there is the choice of gauge- invariant vari- 
ables: the work of Bardeen (1980) made clear that there is no unique choice. Second, 
the use of dimensionless variables invariably leads to physical insight via the choice 
of a suitable normalizing factor or factors. Third, there is the choice of how to 
formulate the governing equations: Einstein's field equations, the Ricci and Bianchi 
identities, the conservation equation for stress-energy, or other matter equations. In 
this paper we systematically consider these three choices, working for the moment 
within the framework of linear perturbation theory. 

Our first goal in this paper is to present an efficient way of constructing di- 
mensionless gauge invariants associated with the metric tensor, the stress-energy 
tensor, or other structures that may be introduced, and of determining their inter- 
relationships. We use the method of Nakamura (2003), adapted as in UW to create 
dimensionless gauge invariants. 

Our second goal is to give a unified treatment of the various ways of writing the 
governing equations in dimensionless form using gauge-invariant variables within 
the framework of the metric-based approach to cosmological perturbations!! In UW 
we gave the linearized Einstein equations in two forms, which we referred to as the 
Poisson form, associated with the work of Bardeen (1980), and the uniform curvature 
form, associated with the work of Kodama and Sasaki (1984). In the present paper 
we derive the linearized conservation equations for the stress-energy tensor and by 
expressing them in terms of suitable gauge invariants, give an alternative description 
of the dynamics of scalar perturbations as a system of two first order (in time) 
partial differential equations. We also include the case where the source has multiple 
components. In addition, by using the inter-relationships between the different 
gauge invariants, we are able to give a unified description of the various " conserved 
quantities" that are associated with long wavelength scalar perturbations. 

Our third goal is to elucidate the connection between the metric-based approach 
and the so-called 1+3 gauge-invariant approach to cosmological perturbation^, 
which was developed with the goal of circumventing the gauge difficulties associ- 

^We follow the nomenclature of Wainwright and Ellis (1997): a Friedmann-Lemaitre (FL) cos- 
mology is a Robertson- Walker (RW) geometry that satisfies Einstein's field equations. 

^By this we mean the standard approach to cosmological perturbations in which one formulates 
the governing equations in terms of gauge-invariant variables associated with the perturbed metric 
tensor and the perturbed stress-energy tensor, using local coordinates. 

^The 1-1-3 gauge-invariant approach uses variables that are apriori gauge-invariant at first 
order due to the Stewart- Walker lemma and are defined using the 1+3 covariant description of 
GR, which is based on a preferred timelike congruence. This approach is growing in popularity. 
For a recent treatment in depth we refer to Tsagas, Challinor and Maartens (2008). 
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ated with scalar perturbations (Hawking (1966), Ellis and Bruni (1989) and Ellis, 
Bruni and Hwang (1990)). The 1 + 3 approach is formulated independently of the 
metric-based approach^ and indeed there is a significant gap between the two ap- 
proaches. In the metric-based approach it is customary to expand the metric and 
other basic variables in terms of a power series in a perturbation parameter as in 
UW, since this clarifies the linearization procedure and permits one to extend the 
analysis to higher order perturbations. In this respect the metric-based approach 
is analogous to standard elementary perturbation procedures in physics and engi- 
neering. On the other hand, the 1-1-3 approach is not formulated as a conventional 
perturbation procedure, and relies instead on deriving exact evolution equations 
which are then linearized by dropping products of first order terms. In this paper 
we will reformulate the 1 + 3 approach so as to bridge the above-mentioned gap. 

The plan of the paper is as follows. In section [2] we give the metric and stress- 
energy gauge invariants and specify the four gauge choices and the two normaliza- 
tions that we will use. In section [3] we discuss the equations for scalar perturbations 
that arise from the linearization of the conservation law for the stress-energy ten- 
sor, and their relation with the linearized Einstein equations. The details of the 
derivation, which makes use of the Replacement Principle in Appendix |Al are given 
in Appendix |Bl In section H] we give a concise derivation of the so-called conserved 
quantities in gauge-invariant form. In section [5] we introduce the basic variables in 
the 1 + 3 gauge-invariant approach to scalar perturbations and derive the governing 
equations, which we then relate to the corresponding equations in the metric-based 
approach. The details are relegated to Appendix O Section E] contains a brief 
summary and discussion. 



2 Gauge invariants and gauge fields 

We begin by describing a dimensionless version of Nakamura's method for construct- 
ing gauge invariants (see Nakamura (2007), equations (2.19), (2.23) and (2.26), and 
UW, section 2.1, for a brief introduction). Consider a family of tensor fields A{e) and 
a background scalar A having dimension length such that X'^A(e) is dimensionles^. 
The change induced in the first order perturbation ^^^A by a gauge transformation 
generated by a dimensionless vector field ^ on the background can be expressed 
using the Lie derivative £: 

A^^U = (1) 

(see, for example, Bruni et al (1997), equation (1.2)). Let X be a dimensionless 
vector field that satisfies 

^X°- = C- (2) 
It follows that the dimensionless object defined by 

A[X]:= A" ((^U - , (3) 



''in the 1+3 approach the metric tensor is not used as a dynamical variable and local coordinates 
are not introduced, in contrast to the metric-based approach. 

^See UW, footnote 9, for a discussion and references about allocation of dimensions. 
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is gauge- invariant. We say that A[X] is the gauge invariant associated with ^^'^A by 
X -compensation. Since a choice of X yields a set of gauge-invariant variables that 
are associated with a specific fully fixed gauge we refer to X as the gauge field. In 
this paper we will use two choices for the normalization factor A: if A is a geometric 
quantity, we will use X = a, where a is the background scale factor, while if A is a 
matter quantity we will use X = Ai, where Ai is definec|§ by f|T5|) . In the latter case 
we will denote A[X] by A[X]. 



2.1 Metric gauge invariants 

Given a 1-parameter family of metrics gab{^)i where e is a perturbation parameter 
and 5'af)(0) is a Robertson- Walker (RW) metric, we define a dimensionless conformal 
metric (jabi,^) according to 

5'afe(e) = a^gabie), (4) 
where a is the scale factor of the RW metric. We expand cjabi^) in powers of e: 

9ab{e) = ^""^gab + e ^^^Qab + • • • , 



and label the unperturbed metric and (linear) metric perturbation according to 



(0)- - few J! (1)- 9gab 

lab ■■= ^ 'gab = 9ab{0), fab ■= ^ ' 9ab 



de 



(5) 



e=0 



In order to construct a gauge field X that satisfies ([2]), using only the metric, we 
need to decompose the metric perturbation fab into scalar, vector and tensor modes. 
Relative to a local coordinate syste nfl we introduce the notation 

/oo = -2v?, (6a) 
fo^ = B,B + Bi, (6b) 
= -2V^7,, + 2D,D,-C + 2D(iCj) + 2Q,-, (6c) 

where the vectors Bi and Ci and the tensor Cij satisfy 

WBi = 0, = 0, C\ = 0, U'dj = 0, 

where Dj is the spatial covariant derivative associated with jij. We can satisfy the 
spatial part AX* = C,"^ of the requirement ([2]) by choosing 

X, = BiC + Ci (7) 

(UW, section 2.2), which we will take to be our default choice for Xj. With this 
choice, the components of the gauge invariant fabfX] associated with the metric 
perturbation fab by X- compensation, are given by (UW, equations (21), (23) and 
(25)) 

foo[X] = -2$[X], (8a) 
foaX] = D,B[X] + B,, (8b) 
f,,-[X] = -2^[X]7i, +2Q,-. (8c) 



^This choice is motivated in section [2?2l 
^See UW equations (10) and (11). 
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where 

^X]:=ip-{drj + H)X'^, ^[X]:=^ + HX^, B[X] := B - 8^,0 + X° , (9a) 

Bj := Bi — driCi, Cij := Cij. (9b) 
In equation fl9al) . 7i is the dimensionless Hubble scalar, defined bj|§ 

n := -. (10) 
a 

The quantities Bj and Cij, which describe the vector mode and tensor mode 
of the perturbation respectively, are intrinsic metric gauge invariant^ that are in- 
dependent of the gauge field X. In contrast the gauge invariants $[-^], "^[X] and 
B[X], which describe the scalar mode, depend on the choice of X° but not on the 
choice of the spatial gauge field X\ Of course, if we leave X* arbitrary then fafe[X] 
contains additional terms and its components are given by 

foo[X] = -2<l>[X], (11a) 

fo^[X] = D,B[X] + B, + d^Zi[X] , (lib) 

fij[X] = -2^X]^ij + 2C,j + 2D(,Z,)[X] . (11c) 

where 

Zi[X] := BiC + Ci - Xi. (12) 
Our default choice ([7]) for X* corresponds to ZJX] = 0. 



2.2 Stress-energy gauge invariants 

Consider a stress-energy tensor T%{e) that obeys the background symmetries, i.e., 
it is spatially homogeneous and isotropic: 

D/°>r'> = 0, (°>r° = (°>r*o = o, ^^>r) = 1 6) (°>r\. (13) 

We assume that T%{e) satisfies the conservation law '^VfeTa(e) = 0, which at zeroth 
order yields 

(oy = -m + ' (14a) 

where 

(o)p = _(o>rOo, = |(°>r\. (14b) 

When constructing dimensionless gauge invariants it is necessary to choose a nor- 
malization factor. In Newtonian theory the dimensionless quantity 5p/ p = ^^^p/^^^p, 
where p is the mass density, is used to describe structure formation. By analogy 
the same quantity is usually used in GR, but with p being the mass-energy density 
instead. We propose that in GR a more natural normalization factor is the inertial 

^Here and elsewhere we denote the derivative of a function /(?]) that depends only on rj by 

^Intrinsic gauge invariants were defined in UW, section 2.1, as gauge invariants constructed 
solely from a single tensor, in contrast to hybrid gauge invariants that are constructed from several 
tensors. 



2 GAUGE INVARIANTS AND GAUGE FIELDS 



6 



mass-energy density ^'^^p + ^^^p, since this is the quantity that appears instead of 
in the relativistic energy-momentum conservation equations. As we will show in this 
paper normalizing with ^'^^p + leads to a simpler description of scalar density per- 
turbations when using matter variables. We shall refer to this type of normalization 
as inertia! mass-density normalization or more briefly, as A4 -normalization. 

In order to implement the above idea we assume that the inertial mass-density 
(o)p _|_ (o)p jj^ ( fi4l) is positive, and introduce a normalization factor with dimension 
length, defined by 

M := {%+^yy'^\ (15) 
As in UW we introduce the notation 

AT:=a\% + ^y), C|:=|^, (16) 

and in analogy with U\V0 we define the following intrinsic gauge invariants associ- 
ated with the stress-energy tensoi^ll, using -normalization: 

^^>r) (17a) 
-M^ (Di(i>r°o + 3-H(^>r°) , (17b) 
M'iCp'>r'o + l^'^\), (17c) 



T 
where 



i^yt) ■= (i>r^. - lsY'^\. (18) 

We also introduce the following gauge invariants by X-compensatior0 using Ai- 
normalization: 

T%[X] := ((^>r°o) - 3HX° (19a) 
T° [X] := (i>r° + DiX°. (19b) 

It follows from ffTTbl) and ([19]) that 

T. = _ (D,n[X] + [X]) . (20) 

Note that stress-energy gauge invariants depend only on the choice of not on 
X\ 

In analogy with UW (see equations (50)) we decompose the matter gauge invari- 
ants Tj, Tj, T, T''j[X] and T'q into scalar, vector, and tensor modes and label them 
as follows: 

T^- = D'^n + 2y'=D(,,II,-) + n^-, (21a) 
T, = Di© + 3i, (21b) 
T = f , (21c) 
T°[X] = DiV[X]+Vi, (21d) 
n[X] = -D[X], (21e) 



^"Replace A by T in equations (39) and (40) in UW. 

"'^"'^The a-normalized gauge invariants T% in UW are related to the corresponding A^-normalized 
gauge invariants T% via T% = AtT%, where At = a^i^% + = (a/M)^. 

i^Replace A by T in equations (38a) and (38b) in UW and multiply by 1/At = {M/af. These 
equations also arise from ([3)) with A replaced by T and X = A4. 
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where 

D^n, = 0, fl\ = 0, Bill) = 0, B% = 0, D^Vi = 0, (21f) 

and 

:= D(,D,) - |7,,D2^ := D*D,. (21g) 

It follows from (EO]) and (EID that 

D = D[X] - 3HV[X], Di = -3nYi. (22) 
2.3 Standard choices of the gauge field 

In order to eliminate the gauge freedom in the scalar mode, thereby determining the 
perturbed metric uniquely, we have to fully specify the gauge field X. We fix the 
spatial part X* of the gauge field ab initio as in equation (I7l), leaving the temporal 
part X*^ to be specified. We observe that X° appears linearly and algebraically in 
the definitions of the gauge invariants: 

nX], B[X], B[X], W[X]. (23) 

Note that ^[X] and B[X] are defined by (Ei]), while ©[X] and V[X] are given b}@ 

©[X] = -M'^^^^T^o + 3HX°, V[X] = A^2D-2j3i(i)yO ^ ^o^ ^24) 

as follows from (fT9l) . (I21dl) . (12 lei) and (I21f[) . We can thus determine X° uniquely by 
requiring that one of these four variables be zero. These choices in fact correspond 
to four of the commonly used gauges in cosmological perturbation theory!^ 



(i) Poisson gauge: 

(ii) Uniform curvature gauge: 

(iii) Total matter gauge: 

(iv) Uniform density gauge: 



B[Xp] = 0. (25a) 

^[Xc] = 0. (25b) 

V[Xv] = 0. (25c) 

D[XJ = 0. (25d) 



-'^^In deriving the expression for V[X], we assume, as in UW, that the inverse operator of exists. 
In terms of the (1 + 3)-decomposition of the stress-energy tensor, we can write Y[X] = v + Q + X", 
as foUows from (fT24|) and (fT2Bl) . 

^^These gauge choices and others are discussed, for example, by Kodama and Sasaki (1984), 
Hwang (1991), Hwang and Noh (1999) and Mahk and Wands (2009). In contrast to our approach 
which emphasizes relations between gauge invariants, they define gauge invariants in terms of 
gauge- variant quantities. 
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Determining in this way does in fact satisfy condition ([2]), AX° = This 
has been verified for Xp and Xc in UW (see equation (26)). For the other two cases, 
we need the transformation laws: 

A(i)T°o = -SnM'^f, A(i)T° = -Al-2D,e°, (26) 

which are a consequence of (HdH Condition ([2]) now follows immediately from 



(El, (EScD, (l25d| and ([26]). 

In practice, we will not use the explicit expressions for X° that are defined 
implicitly by equations (125|) . Instead, in order to be able to relate gauge invariants 
associated with different choices of gauge field we introduce a set of transition rules. 
Let X° be a specific choice of the temporal gauge field and let X° be an arbitrary 
choice. The difference 

Z°[X] :=X°-X°, (27) 
is gauge-invariant on account of ([2]). The desired transition rules are as follows: 

^ [X] = ^ [X,] - [X] , B [X] = B [X,] - Z° [X] , (28a) 

©[X] =©[X.] -S-HZ^iX], V[X] = V[X.] - Z°[X]. (28b) 

Equations fl28al) follow immediately from ( l9a|) . while equations (128bp are a conse- 
quence of ([ISD, fl2Tdl) and fl2TiD . 

By inspection of (BSj) we see that the following linear combinations of the vari- 
ables (1231) are independent of the choice of X: 



[^,B] 
[V,B] 



^[X]-HB[X], [^,V] 
V[X]-B[X], [D,V] 
D[X]-3HB[X], [D,^] 



^[X]-HV[X], (29a) 
D[X]-3HV[X], (29b) 
©[X]-3^[X]. (29c) 



We can thus substitute two different choices of X° into any of the X-independent 
expressions in (12^ and equate the results, thereby relating different gauge invariants. 
For example, if we first choose X° = X° in [D, \E'] and then keep X° arbitrary as 
the second choice, we obtain 

©[Xc] =©[X] -3^[X], (30) 
on account of fl25bl) . If we set X° = X° in this equation it follows that 

D[Xe] = -3^[X,], (31) 

on account of fl25dl) . 

The gauge invariant $ [X] is on a different footing from the gauge invariants (123|) 
since it depends on the derivative of X° through equation (pal) . Thus requiring 
$[X] = does not determine X° uniquely and hence does not lead to a fully 
fixed gaug^. Nevertheless, $[X] does have a well-defined transition rule analogous 
to (12S]) . namely 

$ [X] = $ [X.] + (9, + H) Z° [X] , (32) 



""^^Note that the formal similarity between ([T]) and (jS]) enables one to obtain directly from ((T 
without any calculation. 



16 



As a consequence the synchronous gauge contains residual freedom. 
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as follows from fl9ap . By comparing fl32p with f l28p one can construct X- independent 
linear combinations of and the variables in ( l23ll . analogous to ( |29l) . For exam- 
ple, 

[$, V] := + {dr, + n)Y[X], (33) 

is independent of X. There are three other expressions linking $[X] with ^E'[X], B[X] 
and 3[X] that can be written if needed. If we set X = X^ in (15^ and use f l25cp we 
obtain 

$[X.] = + (9, + ?/)V[X], (34) 

which we will use later. 

To conclude this section we note that the gauge invariants and X-independent 
combinations that we have introduced do not exhaust all possibilities, but do serve to 
illustrate an efficient way of defining gauge invariants and determining their inter- 
relationships, which constitutes one of the main results of this paper. A further 
example arises in appendix IC.2t where we make use of another gauge invariant, 
namely the linear perturbation of the Hubble scalar of a timelike reference congru- 
ence, denoted by H[X]. In working with this gauge invariant we find it necessary 
to introduce an X-independent combination involving three gauge invariants (see 
equation f ll39p ). 

Notation 

Our general notation for dimensionless gauge invariants is exemplified by ^^[Xc] and 
V[Xp], i.e. a capital letter, or a bold face letter, or a special font, e.g. V, which 
denotes inertial mass-density normalization, replaces the symbol for an associated 
gauge- variant variable, with the choice of the gauge vector field indicated by a 
subscript on the symbol X. For convenience we will often simplify the notation by 
setting ^ffX,] = etc.. For some of the commonly used gauge invariants we will 
use unsubscripted symbols: 

$:=$[Xp], ^:=^[Xp], V:=V[Xp], (35a) 

A:=<l>[Xc], B:=B[Xc], © := ©[X,]. (35b) 

3 Structure of the linearized governing equations 

In this section we give different forms for the governing equations for scalar pertur- 
bations, first using the metric gauge invariants as basic variables, and then using 
the stress-energy gauge invariants. 

3.1 Linearized Einstein field equations 

As shown in UW there are two natural choices of intrinsic metric gauge invariants 
when formulating the linearized Einstein equations for scalar perturbations, the 
uniform curvature gauge invariants and the Poisson gauge invariants. As in UW we 
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introduce the geometric background scalars Ac and Cq, with Cq defined in terms of 
the derivative of Ac- 

Ag:=2{-W + H' + K), = -{I + ^CDUAg, (36) 

(see UW, equation (42)). With At and C"^ defined by f|T6l) . the background Einstein 
equations imply that Ac = At and Cq = C"^. We denote their common values by 
A and C^: 

A = Ag = At, C^ = Cl = C|. (37) 



The uniform curvature formulation 

The governing equations in the uniform curvature formulation arj^: 

£bB + a = -n = -^^n (38a) 

H[CaA + C^D^B) = \T + iD^n = At{\T + iD^n), (38b) 

H(D2 + 3ir) B = -iA = -i^T©, (38c) 

HA + {\Ag-K)B= -\V =-\AtN, (38d) 

where the first order differential operators Ca and Cb are defined by 

CA-=d^ + nB, Cb ■■= dr, + 2n, (39a) 

with 



B:=^ + l + 3C^. (39b) 



For future reference we note that 



as follows from (136 



The Poisson formulation 

The governing equations in the Poisson formulation are (UW, equations (54)): 

^ - $ = n = ^tH, (40a) 

(i:-c:^D2)^ = |r+ (iD2 + H£A)n = ^T(|f + (|D2 + Ha^ + 2H')n), (40b) 

(D^ + 3fC)^ = iA = \At^, (40c) 

8^^! + %^= -\V =-l^yV, (40d) 

where the second order differential operator L is defined by 

£(.) := UCaCb (^) , (41) 

^^See UW, equation (52). We give the matter terms on the right hand side in two forms: using 
a-normaUzation as in UW, and A^-normaUzation as introduced in the present paper. 
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or equivalently 

£ = 9j + 3 (1 + C^) ndr, + n^B - (1 + 3C^a)K, (42) 

(UW, equation (56)). For a discussion of these two systems of governing equations, 
and the ways in which they differ, we refer to UW, section 3.2, following equation 
(56). 

3.2 Linearized conservation equations, without Einstein's 
field equations 

As shown in Appendix [Bl linearizing the conservation law VtT^a = leads to the 
following gauge-invariant equations: 

dr^{B[X] - 3^[X]) + D2(V[X] - B[X]) = -3HT, (43a) 

{dr, + n)Y[X] + + = -f - H, (43b) 

where 

S := |(D2 + 3ir)II, (43c) 

and B[X],Y[X], T and 11 are defined by equations (fTTD and (|2T|). 

These equations are valid for any choice of temporal gauge field X°. Referring 
to (l29l) and (133|) we recognize the three groups of terms on the left side as the X- 
independent expressions [ID, [V, B] and [$,V] in (^^. We choose X° = X° and 
X° = X° in the first two, and X° = X° in the third, and use Equations 
then assume the concise form 

9^Dc + D^V = -3Hf , (44a) 
$v + = -f - S, (44b) 

where = D[Xc], V = V[Xp], D = D[X^] and = $[X^], in accordance with our 
convention for labeling gauge invariants. In certain circumstances, the first equation 
can be interpreted as a conservation law for ©c, as will be discussed in section HI 
The second equatioij^ shows that for a barotropic perfect fluid $v is proportional 
to D, and is in fact zero for dust. 

One would hke to use equations to obtain a system of evolution equations 
for the stress-energy gauge invariants D[X] and V[X], for some choice of the gauge 
field X. This is not possible due to the presence of the metric variables 9^\I^[X] 
and $[X]. However, in the case that the stress-energy tensor is the total stress- 
energy tensor one can use the linearized Einstein equations to eliminate these terms 
and achieve the desired goal, as we will show in section 13. 3[ On the other hand 
equations fH5]) are valid for each (non-interacting) individual stress-energy tensor of 
a multi-component source, and as such they form a convenient starting point for 
the derivation of a simple system of governing equations for scalar perturbations of 
such a source. We will derive these equations in section 13. 4[ 



This equation corresponds to equation (5.20) in Bardeen (1980). 
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3.3 Linearized conservation equations in conjunction 
with Einstein's field equations 

We now use the results of sections 13.11 and 13.21 to derive a system of governing 
equations in the form of a first order (in time) system of partial differential equations 
with the stress-energy gauge invariants as basic variables. 

Choose X° = X° in (03]), and eliminate Dp using D = Dp - 3'HV, which is 
obtained from the X-independent invariant [D, V] in ( l29l) . The resulting equations 
are 

d^{I]) + 3HY) -3d^^ + 'D^Y = -3HT, (45a) 
{d^ + n)Y + $ + C|D = -r - S, (45b) 



using the notation ( l35l) . The combination f l45ap — 3'H( ]45bl) can be rearranged to 
read 

(dr, - mCl)B + (D^ + 3K)V - 3(9^^ + 'H^ + \AgY) = SUE, (46) 

where Ac is given by (!36|) . If the stress-energy tensor is the total stress-energy 
tensor, and if we impose Einstein's field equation (140dp and the background field 
equation Ac = At, then the above equation simplifies to 

{dr, - 3HC|)D + (D^ + 3K)Y = SHE. (47) 

Equations ( 1471) and (145bl) form a coupled first order system of evolution equations 
for D and V. However due to the appearance of the metric potential $ the system 
is not closed. We can remedy this deficiency by applying the operator + 3K 
to fl45b|) and using 

Z := (D^ + 3K)Y, (48) 

as a new variable to replace V in the system. On using the Einstein equations f l40ap 
and ( 140d[) . which yield 

(D^ + 3K)^ = (D^ + 3K){-^ - Arfi) = lA{B - 3S), (49) 

equations (HZD, (iHD and (D^ + 3K){<^M>) result in 

{dr, - 3C^n)B + Z = 3nE, (50a) 
[dr, + n)Z+{lA + C^D^ + 3i^))D = -(D' + 3K){T + E) + f^S, (50b) 

where A = Ac = At- For the reader's convenience we note that the variables in 
these equations are defined by equations f l21bp . fl21cl) . fl43cl) and f H8l) . 



Equations fl50|) constitute one of the main results of this paper. They form a cou- 
pled system of first order (in time) partial differential equations for (D, Z), assuming 
that the stress-energy terms T and S are given. They determine the behaviour of the 
scalar mode of linear perturbations of an FL cosmology with arbitrary stress-energy 
content. The structure of this system is similar to the structure of the system of 
evolution equations fl38ap and (138bp for the uniform curvature metric gauge invari- 
ants A and B, and can be derived from them as follows. First use (138dp to express 



A in terms of V. Then apply the operator + 3K to both equations and use fl38cp 
to express (D^ + 3K)B in terms of At^i, after which some obvious manipulations 
lead to equations f l50|) . 
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The evolution equation for © 

By eliminating Z from equations (1501) one can obtain a second order evolution equa- 
tion for the gauge-invariant density perturbation D. We apply the operator drj+Ti to 
the first of equations (l50ll and use the second equation to eliminate Z. The resulting 
equation can be written in the form 

{Cv - C2d2)D = 2(D2 + 3K) (|r + (iD^ + -^q^ + 2'H')n) , (51a) 

where 

:= 9j + (1 - 3C')H9^ + (1 - 3C')H' - (1 + 3C^){n'' + K)- 3{C''yn. (51b) 

Equation f l51al) can also be derived from the governing equations ( l40l) in Poisson 
form. We apply + 3K to fliObl) and use fliOcj) to relate (D^ + S/sT)^ to D. By 
comparing the resulting evolution equation with (15 lap , we can conclude that the 
operator Ct> is related to the operator C according to 

C{A*) = AUi*). (52) 

This result can also be verified by direct calculation. 

Equation fl51ap is a second order linear partial differential equation for D, assum- 
ing that the stress-energy terms T and 11 are given. It differs from other related equa- 
tions in the literature, for example, Ellis et al (1990) (see their equation (48), with 
the coefficients given by equations (19) and (20)), and Hwang and Noh (1999) (see 
their equation (45)), since we have defined D by normalizing with Ad^ = ^^^p + ^^^p, 
while the usual practice is to use If the source is a perfect fiuid with a linear 
equation of state and a cosmological constant, i.e. 

P = Pm + A, p = Pra- A, Pm = WPm, (53) 

then the right side of (15 lap is zero, and one can use Einstein's equations in the 
background model (UW equations (41)) to write the expression fl51bp in the formic 

Cv = + (1 - 3w)ndr, - |[(1 + 3w)(l - w)pm + AwA]a^. (54) 

In this case (15 lap is compatible with those equations cited above. 

3.4 Governing equations for a multi-component source 

The governing equations for a perturbed PL cosmology with a multi-component 
source were first derived by Kodama and Sasaki (1984), and subsequently considered 
by various authors including Hwang (1991), Dunsby et al (1992) and Durrer (2008). 
These authors, as is customary, use normalization with ^^^p when defining the density 
perturbation. We have found that the derivation and the form of the governing 
equations is significantly simpler if one uses Al-normalization, as introduced in 
section 12.21 In this section we thus give a brief derivation of the relevant equations. 



Since w = constant we have — w. 



3 STRUCTURE OF THE LINEARIZED GOVERNING EQUATIONS 



14 



We consider a multi-component source with n separate stress-energy tensors 
denoted by aT%, with A = l,...,n, which sum to form the total stress-energy 
tensor: 



T% = Y,aT%. (55) 



A=l 



For simplicity we assume that the individual components are non-interacting. As 
shown in Appendix |Bl linearizing the conservation equation ^ aT\ = for an 



arbitrary component labeled A leads to the following equation^ 

dr,(PA[X] - 3^[X]) + B\Wa[X] - B[X]) = -SnVA, (56a) 

(dr, + n)YA[X] + $[X] + CIBa = -Ta - Ha, (56b) 

where 

S^:=|(D2 + 3K)n^. (57) 

We assume that the gauge field X does not depend on the labeling index. A quantity 
¥a associated with aT% that satisfies 



J2I3a^a = ^, (5^ 



A=l 

where F is the corresponding quantity associated with T%, will be called additive. 
Here the coefficients Ba are defined by 

Aa " 
Ba '■= — j-, and satisfy ^^^Ba = 1- (59) 

A=l 

We note that Ba[X], Ya[X], Ba, ^a, f a + C^Ba[X] and Cl are additiv(|l]. 

In order to obtain a closed system of evolution equations we introduce the " dif- 
ference variables" 

Bab := Ba[X] - Bb[X], Yab := Ya[X] - Yb[X]. (60) 

It follows from ffTOal) . ffT9bl) . fl2Tdl) and fl2lel) that Bab and Yab are X-independent. 
In order to obtain evolution equations for Bab and Yab we form the difference of 
two copies of equations (l56ll . labeled A and B. In this way we obtain the following 
equation^: 

dr,BAB + D^Vab = -S-HTab, (61a) 
(a, + H)Yab + Kab + {Cl - Cl)B = -Tab - ^ab, (61b) 



these equations denotes the value of for the component labeled A. We have dropped 
the subscript T for convenience. 

^^One can check the consistency of (|56l) by showing that y^^_-^ ,64 ((56)) = p3|) . Note that 
]S])A='BA[X]-3nYA[X]. 

^^The only calculation involves showing that C^Da — C^Ds = Kab + {C\ — C^)©, which follows 
by writing = D + X]c=i ^c{^A — ©c), and a similar expression for Ds. 
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where 

n 

Kab := Yl [Cli^Ac - mNAc) - CliJ^BC - 3HVbc)) , (61c) 

C=l 

Tab := Ta — F^, S^b := 2^ — S^, (61d) 

Equations ( 16T|) do not form a closed system for D^b and V^b, however, due to the 
appearance of the total intrinsic matter gauge invariant D. The evolution equa- 
tion (15 lap for D contains F and S as source term@, which have to be expressed in 
terms of f a, and D^b- The term f is not additive, whereas f + C^D[X] is, i.e. 

n 

f + C2d[X] = Y,^A{fA + Cl'BA[X]). (62) 

A=l 

It follows that f can be expressecQ as a sum involving f ^ and Dab: 

n n 

f = ^ BaTa + \ ^aBb{CI - CD^AB. (63a) 

A=l A,B=1 

On the other hand, the term S is additive: 

n 

E = Y13a^a. (63b) 

A=l 

In conclusion, equations ( l6Ti) and the evolution equation f l51al) for D, in conjunc- 
tion with equations (163|) . form a closed system for ]D)ab,^ab and D that describes 
the scalar mode of a perturbed PL cosmology with a multi-component sourcj^. We 
note that all gauge invariants in these equations are X-independent. 



3.5 Gauge-fixing versus gauge- invar iance 

To conclude this section we comment briefly on the two points of view as regards 
formulating the governing equations for scalar perturbations. In the gauge-fixing 
approach one chooses a gauge ab initio, in which case all the gauge invariants that 
appear in the governing equations are associated with this particular gauge. Recent 
references that use this traditional approach are Mukhanov (2005) and Weinberg 
(2008). In contrast, one can work with a variety of gauge invariants in the spirit 
of Bardeen (1980), in which case one has the flexibility to use gauge invariants as- 
sociated with different gauges in formulating the governing equations. A notable 

^■^Note that the term (D^ + 3K)Il on the right side of (|51ap can be replaced by S on account 
of (143^ . 

24Substitute Ba[X] = B[X] + Y.c=i ^c^ac in (US) and use J2a=i ^aCa = 
^^For the reader's convenience we give the equation numbers in the previously mentioned ref- 
erences that correspond to our equations (|6T|) : Kodama and Sasaki (1984), (5.59)-(5.60), Hwang 
(1991), (37)-(38), Bruni et al (1992b), (91)-(92), Dunsby et al (1992), (86)-(87), and Durrer (2008), 
(2.136)-(2.137). 
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example is the generalized Poisson equation f l40cp . which relates the Bardeen po- 
tential \1/ = ^I/[-^p] to the matter density gauge invariant D = D[Xv]. A variety of 
other examples occur in this paper, including the system of equations f l50|) for ©v 
and Zp, equations f l38cp and (I38dp in the uniform curvature formulation, and the 
conservation law (l66l) . 



4 Conserved quantities 

Two "conserved quantities" that are associated with scalar perturbations of FL 
have been defined in the literature. These quantities, often denoted by (, satisfy an 
evolution equation of the form 

a,C. = D2C. + ?^f, (64) 

where C, is an expression involving the primary gauge invariants such as ^ or V 
and the background variables. This equation is referred to as a "conservation law", 
since if spatial derivatives are negligible ("perturbations outside the horizon") and 
if r is zero or negligible in some epoch, then flM]) is approximated by 9^C» = 0? '^•C- 
C, is approximately constant in time during that epoch. 

Two of the evolution equations that we have presented in section [3l namely 
equations f l44ap and (]38bp . can be written in the form ( |6^ simply by multiplying 
by a suitable factor. The conserved quantity can then be identified by inspection. 
We consider each in turn. 

The conserved quantity (p 

The evolution equation f l44ap for Dc has the form of a conservation law (!64l) . We mul- 
tiply by the numerical factor — | to agree with current convention, and comparison 
with leads to the following definition of the first conserved quantity: 

Cp := -l^c (65) 

Equation fl44al) assumes the form 

d^Cp = iD^v + nr. (66) 

Our motivation for using the notation (p is that on account of equation (13T!) we have 

Cp = %■ (67) 

This conservation law was apparently first given in a form closely related to the 
above by Wands et al (2000) o "who emphasized that it depends only on the con- 
servation equation for the stress-energy tensor, i.e. it is independent of Einstein's 
equations. They denoted our Cp by C and because of ( E7|) they referred to it as "the 

^^See their equations (8) and (9). Their evolution equation (18) corresponds to our equation ([TO)) , 
although their equation is not in a manifestly gauge- invariant form, and uses clock time rather than 
conformal time. See also Malik and Wands (2009), equations (7.61), (7.62) and (8.35), which are 
somewhat closer in form to our equations. 
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curvature perturbation on uniform density surfaces." The quantity (p has its origins 
in the paper Bardeen, Steinhardt and Turner (1983), and was further studied from 
a different point of view by Brandenberger and Khan (1984)0. 

A major step in understanding the significance of (p was taken by Langlois and 
Vernizzi (2005). Motivated by the 1 + 3 covariant approach, they showed that 
this quantity could be obtained as the hnearization of an exact nonhnear evolution 
equation for a certain covector. This approach enabled them to extend the definition 
of (p to second-order (and higher order) perturbations. We refer to their equations 
(20) and (25) for the general situation and note that their equations (41) and (42) 
correspond to our equations (IB7|) and (^E^. 



The conserved quantity 

The second conservation equation arises from the evolution equation fl38b|) for A = 
$c- Using fl39cp we can write the differential operator Ca in the form: 



/2 



HCa{*) = 4^0, ] . (6^ 



Thus on multiplying ( I38b|) by 21-1/ Ac we can write it in the form of a conservation 
la^ 

d, {^^^ = (^C'j + in^ + nf, (69) 

where we have chosen to replace B by using the relation \E' = — "HB. Compar- 
ing (169|) with fl6^ leads to the following definition of the second conserved quantity: 

Cv := ^$c. (70) 

Equation (1691) . which is equivalent to fl38b|) . is the conservation equation for An 
immediate consequence of this equation is that if the source is pressure-free matter 
plus possibly a cosmological constant, then Cv is constant in time. 

We can derive an alternate expression for the conserved quantity as follows. 
Using the relation = — 7/B and the definition fITO]) . the velocity equation fl38dp 
can be written in the form: 

Cv=fl-^V-^V. (71) 



A J 

We use the X-independent gauge invariant [\1', V] in (12^ to obtain 

= ^ - -HV, (72) 



^^It is not immediately obvious that the expressions given in these papers (equations (2.43) and 
(2.45) in Bardeen et al (1983) and equations (2.11) and (2.12) in Brandenberger and Khan (1984)) 
agree with our expressions. 

^^Throughout the remaining discussion about conserved quantities we use the background Ein- 
stein equations and hence Ac = At = A and Cq = — C^. 
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which when inserted in fl7ip leads to 

2K 

Cv = - -j^. (73) 

The quantity is most commonly used when the background curvature is zero 
(i^' = 0) in which case 

Cvk=o = (74) 

This expression motivates our use of the notation Malik and Wands (2009) use 
the notation 71 for in the case K = (see equation (7.46)), and refer to it as "the 
curvature perturbation in the comoving gauge." This quantity has its origin in the 
paper Bardeen (1980) (see equations (5.19) and (5.21)). 

There is another commonly used expression for ^v, in terms of the Bardeen 
potential \1/ and its time derivative, which we can quickly derive. On solving (I38ap 
for A = the definition (170|) yields 

c. . (c. (I) - n) . (75) 



Ac \ \n, 

where we used "HB = — \I' and where £b is defined by fl39al) . We expand the operator 
and use (15^ to express Ti' in terms of Ac- On specializing to flat FL {i.e. K = 0) 
and setting 11 = we obtain the expression 

^-i-» = * + 3(rfT;)(w*'+*)- 

where w := ^^)j/^'^'>p. Here we have used the background field equation Ac = At, and 
the fact that At = 3(l+w)7{^ when = 0, as follows from UW (see equations (41a) 
and (42a)). The familiar expression fl76l) can be found, for example, in Mukhanov 
et al (1992), equation (5.23), and Mukhanov (2005), equation (7.73). 

We conclude this section by showing that Cp and Cv) despite their different origins, 
are closely related. On account of fl67j) and (!73l) 

2K 2K 
Cp - Cv = - + -^^ = -1© + (77) 

where we have used the X-independent invariant [D, in ( 129|) to obtain the second 
equality. The generalized Poisson equation fl40cp can be used to eliminate D yielding 

Cp - Cv = -I^D'*- (78) 

This equation suggests that if spatial derivatives are negligible in some epoch, then 
Cp ~ Cv in that epoch. 

A coupled system for ((^^, ^) 

The conservation equation (1690 for is one of the linearized Einstein equations in 
uniform curvature form, namely, the evolution equation for A = $c- It is helpful to 
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also write the evolution equation f l38ap for B in terms of Cv, replacing B by \1/ using 
the relation = — T/B. The resulting pair of equations has the following form: 



5^(^-^J -H. = x'AU, (79a) 
9,Cv - yD^ (^vl/^ = nif + fD^n) , (79b) 



where x := a/a* is the dimensionless scale factor, with a* being the scale factor at 
some reference time, and 

The system of equations (1791) is a particularly useful form of the governing equations 
for scalar perturbations of FL. The second equation is the "conservation equation" 
for Cv, while the first equation enables one to express the Bardeen potential \1/ as a 
quadrature, in situations in which is a temporal constant (or can be treated as 
such) and 11 is negligible: 



vl/(r/,xO = ^ (Cv(xO^VrffJ + C_(xO) 



In particular this equation gives the exact solution where the source is pressure-free 
matter and, possibly, a cosmological constant (since then is a temporal constant), 
and the approximate solution in the long wavelength limit, in both cases without 
restriction on the background spatial curvature K. 



5 The 1 + 3 gauge-invariant approach 

In this section we first give a concise derivation of the governing equations for linear 
scalar perturbations in the 1 + 3 approach, combining the formulation of Bruni, 
Dunsby and Ellis (1992a) (hereafter referred to as BDE)owith our overall strategy of 
creating dimensionless quantities by normalizing with Ai and a. We then introduce 
the dependence of the 1 + 3 variables and differential operators on a perturbation 
parameter e, which enables us to relate the variables and governing equations of the 
two approaches in a precise manner. In this way we set the stage for extending the 
1 + 3 approach to second order. 

The 1 + 3 gauge-invariant approach to cosmological perturbations is based on 
choosing a preferred unit timelike vector field and decomposing the stress-energy 
tensor relative to this vector field: 

n = (P + PKuk + p6% + {q'^Uk + u'^q,) + n, (82) 

where 

M> = 0, 7r'^a = 0, UaTl% = Q. (83) 

^^BDE give a comprehensive account of the hnearization of the full system of equations in the 
1 + 3 formalism (Ricci and Bianchi identities, and stress-energy conservation equations ). We are 
concerned only with a limited subset of these equations. 
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One distinguishes between physical and geometrical quantities which are non-zero in 
the background spacetime, namely p,p and the Hubble scalar H associated with m", 
and quantities which are zero in the background, such as the stress-energy quantities 
qa, 71%, the shear of the preferred congruence, the Weyl curvature tensor and the 
spatial gradients of p,p and H orthogonal to u"". The 1 + 3 approach focusses on 
the latter quantities, which are gauge-invariant on account of the Stewart- Walker 
Lemma. 



5.1 Evolution equations 

The spatial gradients of p and H describe the perturbation in a gauge-invariant way 
at the linear level. However, in order to extract the scalar mode of the perturbation 
it is necessary to form scalar quantities. We thus take the spatial divergence of these 
spatial gradients and form the dimensionless spatial Laplacian: 

D := (a' ^^^V^)M^p, Z := ?,{a^ ^''W^)aH, (84) 

where 

V, := K'V,, := g^' ^'^Vj'W,, (85a) 

and 

ha' := 5a' + UaU\ (85b) 

Note that in introducing dimensionless variables we are normalizing the energy den- 
sity p with Ai^ and the geometric quantity H with the background scale factor a. 
Likewise we normalize the geometric operator ('^^V^ with a^. 

The governing equations for the scalar mode take the form of a coupled system 
of first order (in time) partial differential equations for D and Z. These equations 
arise from the energy conservation equation (the exact evolution equation for p) 
and the Raychaudhuri equation (the exact evolution equation for H). To derive 
the governing equations one simply applies the differential operator to the 

linearized versions of the evolution equations for p and H, which are obtained by 
dropping products of first order quantites. 

The linearized evolution equations for D and Z, derived in Appendix [Cl are as 
foUowS 

D' - 3nC^D + Z = 3H(n + t) - D^Q, (86a) 
M + C|(D2 + 3K)\d = + 3K)(T + U + f) + |^(n + T), (86b) 



z' + nz + 



In these equations the dimensionless operators ' and are defined by 

A' := au^'VaA, &A := ^^'^V^A, (87) 



^^Equations equivalent to (1551) have been derived by Woszczyna and Kulak (1989) in the case 
of a barotropic perfect fluid, using the method of Appendix [C] but with different normalization 
factors (see their equations (11) and (18)). Their variables Ae and A9 are related to D and Z 
according to D = a^M'^Ae, Z = a^AO. 
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where A is a scalar. The source terms 11, Q and T are first order dimensionless 
scalars formed by taking the spatial divergence of and 7i% after normalizing with 

Q:=I)^{M\a), n := D,D^(MVfc), f := Q' ~ {3C^ - l)nQ, (88) 
where 

:= a^^^Va. (89) 
The entropy perturbation f is given b\l^ 

f = P- C^D, (90) 

where 

P := I)\M^p). (91) 

We conclude this section by relating our approach to that of BDeI^. The variables 
D and Z differ from those introduced by BDE as regards the normalization of pm 
and H. Specifically, BDE define 

V Pm J 

At the linear level, the factor pm in the denominator can be replaced by '•^^pm- In the 
1 + 3 approach the scale factor is usually defined using the Hubble scalar H of the 
preferred congruence m", according to {u°''Va a)/a = H,ui which case (3)V"a ^ 0. At 
the linear level, however, the factor of a can be taken outside (■^^V", since the term 
(3)Y'"a appears as a product with '•^^V"pm or ^^^'^aH, and hence can be dropped. 
To avoid this complication we have chosen the scale factor a to be the scale factor 
in the background model so that (^)V"a = 0. This choice also facilitates the link 
with the metric-based approach. In view of these remarks, the BDE variables (p2|l 
can be written in the form: 

A = TT^D'pm, Z = SD'if, (93) 



which implies that the BDE variables are related to ours according tco 

A = {l + w)D, aZ = Z, (94) 

since 1/ Ai"^ = {l + w^'^^pm- Our evolution equations ( l86i) for D and Z are equivalent 
to equations (68) and (69) for A and Z in BDE, but are simpler in form due to our 
use of dimensionless variables, in particular our use of -normalization!^. 



3iSee BDE equations (27), (28) and (40). Note that T = 

■^•^BDE do not incorporate the cosmological constant into the stress-energy tensor as we do. In 
making a comparison we have to write p = /9,„ + A, p = p,„ — A, with w — Pm/ Pm not necessarily 
constant. 

■^■^Note that p + p = pm + Pm = (1 + w)pm, and that D^pm = D^p. 
■^^Our variables {Q,II,1l) are related to those used in the above reference according to 
^BZJB = Q, allBDE = n, aFsDE = T. 
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5.2 Relation with the metric-based approach 

Equations f l86|) are closely related to the governing equations in the form fl50|) for 
the variables D and Z, that arise in the metric-based approach. Indeed a formal 
similarity is obvious on inspection. However, in order to relate the two sets of 
equations we have to regard each of the variables D, Z, Q, U and F in ( l86i) as being 
a function of the perturbation parameter e, which can be expanded in a power series 
of the form: 

F(e) = (o)F + e«F + ... . (95) 

Since each variable is zero in the background we have = 0, while ^^^F is the 
linear perturbation of F. We also need to consider the dependence of the differential 
operators on e, which is as follows: 

{A'){e) = au^ieyVaAie), {I)'A){e) = aV{eYVa'V,A{e) (96) 

Assuming that A is a scalar such that ^'^^A = 0, it follows that 

(i)(A') = ^'\-b^A) = (^U, (97) 

as is shown in Appendix [Cl 

If we now differentiate equations fISBl) with respect to e and set e = the resulting 
equations have precisely the same form but with each variable replaced by its linear 
perturbation and each differential operator replaced by the corresponding zeroth 
order operator. We finally have to do a calculation using (!95l) to specifically relate 
the variables in the two sets of equations (!86l) and (l50l) . The details are given in 



Appendix O, where it is shown that: 

D^ro = - m^^'^Q, (98a) 

D^Z = + (D^ + ?>K- |^g)^'^Q (98b) 

D^H = (^ft, (98c) 

D^f = (98d) 



If we now apply the operator to equations ( !50|) and use equations ( !98l) then we 
obtain precisely equations ( l86|) . with each variable by its linear perturbation and 
each differential operator replaced by its zeroth order perturbation. 

6 Discussion 

In this paper we have presented an efficient way of defining dimensionless gauge 
invariants and determining their inter-relationships, which we have applied to give a 
unified account of the various ways of formulating the governing equations for scalar 
perturbations of FL cosmologies. In defining gauge invariants we use our version of 
Nakamura's geometrical method, as described in UW (see section 2.1) which is based 
on specifying a so-called gauge field X, and normalizing so as to obtain dimensionless 
quantities. It turns out that the choice of the spatial part X* of the gauge field does 
not affect the form of the governing equations for linear perturbations given in UW 
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and in the present papered In these papers, however, we find it convenient to fix 
the spatial part X* of the gauge field according to ([7]), which leads to a simple 
form dH]) for the metric gauge invariant faf)[X]. This in turn shortens the calculation 
of the Riemann gauge invariants in UW (see (B.23)) and of the gauge-invariant 
form of the divergence of the stress-energy tensor in the present paper, equation 
f l43p F^ The remaining gauge freedom is then described by the temporal part 
which we specify uniquely by requiring that one of the four basic gauge invariants 
^PfX], B[X], D[X] and Y[X] be zero. This approach eliminates the need to express 

explicitly in terms of gauge-variant variables, and in particular enables one to 
perform a change of gauge without using X°, as in subsection 12.31 Indeed, although 
the gauge field X", which is gauge- variant, plays an important role in establishing 
our formalism, we do not use it in performing calculations, in keeping with our goal 
of working exclusively with gauge invariants. Because it simplifies calculations this 
approach will facilitate the extension to second order perturbations. 

The coupled system of first order partial differential equations for D and Z, given 
by equations fISUl) . plays a central role in this paper. To the best of our knowledge 
they have not been given in the literature. They arise first of all from the linearized 
conservation equations in conjunction with the linearized Einstein equations, but can 
also be derived directly from the latter equations, when they are written in terms 
of the uniform curvature gauge invariants as in (1381) . Further, as shown in section [5] 
these equations are essentially equivalent to the first order governing equations for 
D and Z that arise in the 1-1-3 gauge-invariant approach (see equations (15^ ). 

Our derivation of the expressions for the conserved quantities and deserves 
comment. We have shown that the conservation equations for (p and are simply 
two of the first order governing equations for scalar perturbations when they are 
written in terms of the appropriate gauge invariants: the gauge-invariant expression 
for (p arises directly from the linearized conservation equations for the stress-energy 
tensor, while that for arises from the linearized Einstein equations for the uniform 
curvature metric gauge invariants. Other expressions for the conserved quantities 
are derived using the method for finding inter-relationships between gauge invariants 
given in section 1^31 We mention that Cv is usually introduced by rewriting the second 
order evolution equation ( I40b|) for the Bardeen potential \E' in a first order form, a 



procedure that involves a tedious calculation^] Our use of the uniform curvature 
metric gauge invariants avoids the need for any calculation. 

Our discussion of the 1-1-3 gauge-invariant approach has several novel features. 
An advantage of this approach is that it is coordinate- free, so that calculations re- 
quire only standard operations from differential geometry. This feature has enabled 
us to give a particularly concise derivation of the first order system of governing 



^^This feature of linear perturbations is due to the fact that the components of the perturbed 
Riemann tensor and of the perturbed stress-energy tensor that are used in deriving the governing 
equations are invariant under spatial gauge transformations. Indeed, according to Bardeen (1988), 
since the background 3-space is homogeneous and isotropic the perturbations in all physical quan- 
tities must be invariant under spatial gauge transformations. Whether this property holds for 
nonlinear perturbations requires further investigation. 

■^^If one does not fix X' one has to use the general form (fTTj) of fab[-'^]- However, during the 
calculations the terms involving X* cancel, leading to the same final results. 

^''This calculation is much easier if one uses our factorization property (|4T|) for the operator £. 
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equations for scalar perturbations, as given by equations (l86ll (see Appendix O) 
We have also derived the relation between the variables (©, Z) in the metric-based 
approach and the 1 + 3 variables {D,Z) (see equation fl98|) ). A drawback of the 
1 + 3 approach is that the linearization process is conceptually less clear than in 
the metric-based approach, relying as it does on "dropping products of first order 
terms" . In relating the 1 + 3 approach to the metric-based approach it was necessary 
to regard D, Z and the differential operators as functions of the perturbation pa- 
rameter e and explicitly calculate their dependence on e to linear order. Introducing 
the perturbation parameter clarifies the linearization process and points the way to 
extending the 1 + 3 approach to second order perturbations. 
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A The Replacement Principle 

We define 

4(e) ■.= M''VtT\ie). (99) 

The linear perturbation of la, given in equation flllOp . can be written symbolically 
in the form 

(^)/a = U(A^^«n, (100) 

where Lq is a linear operator. The replacement principle for the divergence of the 
stress-energy tensor states that the gauge invariants associated with ^^^la, ^^^T% and 
fab by X- compensation are related by the same linear operator: 

I„[X] = U(T^,[X],f,4X]), (101) 

for any gauge field X. If the stress-energy tensor is conserved at zero order {i.e. 
-^a(O) = then ^^'^la is a gauge invariant, and the left sides of (llOOp and (11011) are 
equal. 

This result is adapted from Nakamura (2005) (see equations (3.90), (3.91) and 
(3.20)). Use of the Replacement Principle in Appendix [B] makes the transition from 
gauge- variant to gauge-invariant equations particularly easy and transparent. 



B Derivation of the conservation equations 

In this Appendix we give the derivation of the linearized conservation equations in 
the form fH3|) . using the methods developed in UW (see in particular Section 2 and 

^^In the usual derivation one first obtains a system of partial differential equations for the spatial 
gradients of the energy density and the Hubble scalar, and then one takes the spatial divergence 
to obtain partial differential equations for scalars, a more lengthy process. See, for example, BDE. 
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Appendix B). We express the covariant derivative '^'Va of the metric Qabi^) in terms 
of the covariant derivative °Va of the conformal background metric 7^6 as follows: 

'VaA',{e) = °VaA^e(e) + QU^)Ai{e) - QU^)A''4e). (102) 

The object Q%c{e) is written as the sum of two parts: 

Q%{e) = Q%c{e) + QU^), (103) 

where 

Q%c{e) ■■= 25\brc) - r''{e)9U^)rd, with := °V,(lna), (104a) 

Q%c{e) := ir'(e) (°V,^d,(e) - V,^,e(e) + °Vfe^,d(e)) . (104b) 

It follows from fll04al) and fll04bp . in conjunction with ^albc = 0, that at zeroth 
and first order we obtain 

= 25%r,) - r%crd, ^^^Q%c = 0, (105a) 
^''^Q%c = {r'lbc - 7"'/f>c)r,, ^^^Q%c = \ (°Vc fdb - + Ud) . (105b) 

Consider tensors A%{e) such that \'^A%{e) is dimensionless, where A > is a 
background quantity with dimension length. As follows from (11021) . the equation 
= ""VbA a(e) can be written as 

= (°Vb - 2sk) A^^Ue) + 2Q\,{e)\^A%e), (106) 

where 

s,:=°V„(lnA), (107) 
which yields the following zeroth and first order expressions 

= (°V, - 2.,) \'^'^A\ + 2WgVA'^°U^a], (108a) 
= (°V, - 2.,) A^^^U^, + 2WgVA'^'^A^a] + 2(i)g^,[,A2(°U'L]. (108b) 

We now specialize A%{e) to a stress-energy tensor T%{t) that obeys the back- 
ground symmetries, i.e. that satisfies ( fT3l) . and choose the normalizing factor A as 
in equation (fT5|) . ie. A = A^. We also assume that T%{e) satisfies a conservation 
law of the form ^Vf,T ^(e) = 0. Relative to local coordinates we obtain 

Ta = ^5 a, = |-H(1 + C|)5°„. (109) 

On substituting from (11051) the zeroth order expression fll08al) yields equations (fl^ , 
and the temporal and spatial components of the first order expression (]108bp assume 
the following form: 

= d,{M' (^>r°o - ifi) + Bi{M^ (i>r*o) - nM^ {^'>T\ + 3C|(^>r°o) , (llOa) 

= (9, - mcl){M^ (^>r° ) + D,(A^2 ^^>T\) - Hi,, - |D,/oo + nu 

(110b) 
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We simplify these equations by first expressing (i>ri) in terms of 

A^2(i>r^o = -i'' (>i'(^>r° - /o,) , (111) 

and decomposing (^^T^ into its tracefree part and its trace using (fTSI) . We then 
introduce the intrinsic gauge invariants T^i,Ti and T as defined by (ITTI) . expressing 
the trace ^^^\ in terms of T. As a result of these changes equations (II 101) yield: 

= d,{M' - ifi) - B'iM^ - fo^) - 3-HT, (112a) 
= {d^ + n){M^ - ^D,/oo + C^T,, + + D,T. (112b) 

We now apply the Replacement Principle to these equations, which entails perform- 
ing the following replacements: 

/oo^foo[X], /o. ^fojX], /,^. ^f,^.[X], (113a) 

On substituting from ([8]) and (I2T1) and noting that 

D^fo^X] = B^B[X], D*T° [X] = Ti'^Y[X], (114) 

equation (I112ap assumes the form (143 al) . After performing the replacements (I113p 
in ( I112bp we apply the operator D* in order to extract the scalar mode. We then 
substitute from ([8]) and ([21]), noting (I114p and the fact thati^ 

DT, = D^D, B'BjTi = D^D^D^ill = D^S, (115) 

where S is defined by (I43cp . The result is that (]112bp assumes the form D^C = 0. 
Since we are assuming, as in UW, that the inverse operator of exists, we obtain 
C = 0, which is precisely the desired equation (]43bp . 



C Derivation of the 1 + 3 perturbation equations 

C.l Derivation of the evolution equations 
The evolution equation for D 

We begin with the conservation equations for the stress-energy tensor fl52]) . linearized 
by dropping products of first order quantitiej^: 

p' = -3aH{p + p) -B^qa, (116a) 
haq'b = -4:aHqa - D^p - (p + p)aUa - Dfevr (116b) 
■^^The third equality follows from the identity (B.39e) in UW. 

^°See, for example, Wainwright and Ellis (1997), equations (1.48) and (1.49), after multiplying 
by a to change the dot derivative to prime. 
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In these equations the differential operators ' and Da are defined by flHT]) and fl89l) . 
We require the zero order version of flll6ap which we write in the forno 

M^'^y = -3'H, (117) 

which leads to the evolution equation for Ai^: 

{M^y = 3{i + c^)nM\ (118) 

We apply the operator M'^'D'^ to flll6ap and the operator A^^D" to flll6bp and 
then linearize, obtaining 

D' - mC^D + Z = -m{t>''{aua) + P) - D^Q, (119a) 
D"(auJ + P = -(T + n). (119b) 

On substituting flll9bp in flll9ap we obtain the evolution equation fl86ap for D. 

In deriving (I119ap and (I119bp we use the following linearized commutativity 
properties: 

t>\A') = (t>^Ay - (°U'(D'^(aMa)) , (120a) 
where A is any scalar field, and 

D.A; = {BaAtY, (120b) 

where Aa is any covariant vector field. In differentiating products of perturbed 
quantities such as pH, Hqa and (p + p)ua we use the following expansion to linear 
order: 

AB = ^''^AB + ^''^BA - (121) 

where A and B are geometric quantities with background values ^^^A and ^^^B, one 
of which may be zero. 



The evolution equation for Z 

We begin with the linearized Raychaudhuri equation written in the form 

3{aH' + a^H"^) - D"(aMj + ^a^{p + 3p) = 0. (122) 

We use ffml) with A = B = H to write a^H^ = 2H{aH) - H^, where H := a^i/, 
and use ( 1119bp to eliminate Ua- We then apply the operator to f ll22p . After 
using fll20ap and the definitions oi D^Z and P we obtair@ 

Z' + HZ + \AD = -(D^ + 3K)P - (D^ + 3K- |^)(t + H). (123) 

In deriving this equation we have also used fl37|) and fll37p . We finally use fl90|) to 
express P in (11230 in terms of D and F, which gives the evolution equation fl86bp 
for Z. 

^""^For a background scalar ^^^A' is the ordinary derivative with respect to confornial time rj. 
^^In doing calculations such as these one should keep in mind that ''^^Va(a') ^ 0, where ' is 
defined by (|87p. even though we have chosen a to be the background scale factor (i.e. Va(a) = 0. 
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C.2 Relation between the 1 + 3 and the metric-based ap- 
proaches 

Fundamental 4-velocity and energy flow vector 

We begin with the decomposition of the stress-energy tensor with respect to a unit 
timehke vector field u", which is given by (!82l) . The Stewart- Walker lemma implies 
that the linear perturbation ^^\a is a gauge invariant. Since '•^V" = and 
^^^Ua = — a5°a, it follows that ^^•'go = 0, and hence that 

^'>T'o = -%, M'^'yT%=v, + Q,, (124a) 

where 

av, := aQ, := (124b) 
It follows from f lT24|) and (fTOb]) that 

T°[X] = v4X] + Q„ (125a) 

where 

v^X] = i;, + DiX°. (125b) 
We decompose Qj, Vj and Vi according to 

Qi = DiQ + 4, v,[X] = Div[X] + Vi, Vi = D,v + Vi, (126a) 

with 

D*Qi = 0, D^Vi = 0, mi = 0. (126b) 
It now follows from fl2ldl) . fn:25al) and ffT26D that 



V[X] = v[X] + Q, V, = v, + Q,. (127) 

Thus if the preferred timelike vector field is an eigenvector of the stress-energy 
tensor, i.e. if the energy transfer vector is zero, then the stress-energy gauge 
invariants Y[X] and Vj equal the gauge invariants v[X] and Vj associated with m". 
In addition it follows from ( ]125bl) and f ll26p that 



v[X]=t; + X°, ^i=Vi. (128) 
Spatial gradient and Laplacian of a scalar 

We have seen that the 1 + 3 approach to cosmological perturbations is based on 
the spatial gradient and Laplacian of the density p and the Hubble scalar H. We 
now define these quantities for a scalar field of given dimension, using a background 
normalization factor A of dimension length. Let / be a scalar such that A"/ is 
dimensionless, and whose unperturbed value is a function only of rj. We define the 
dimensionless spatial gradient and spatial Laplacian of / according to 

:= 1d,(A"/), F:=&{X-f), (129) 
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using the notation (153]) . 

Our goal is to relate the linear perturbation of Fa and F to the linear perturbation 
of /. Regarding all perturbed quantities as functions of the perturbation parameter 
e, we write 

Fa{e) := K\erVaiX''f{e)), F{e) := aV'(e)/ia'^(e)^V,n(e), (130) 
with /(e) = '•°^/ + e ^^^f + . . . , etc. A straightforward calculation yieldS 

(°)F, = 0, «Fo = 0, «F, = A"(D,«/ + (°)A.)' (131a) 

(0)^ = 0, (^^F = (^^Fi. (131b) 

We note that the background values '•^^Fq and *^°^F are zero due to our assumption 
that *^°^/ = ^^^f{r]). On account of the Stewart- Walker lemma the linear perturbations 
and ^^^F are gauge-invariant. We can write them in a manifestly gauge-invariant 
form by noting that 

A" (D/i)/ = D,f [X] + A" ("^/'v^X], (132) 

where f[X] is the gauge invariant associated with / by X-compensation and Vj[X] 
is given by i^2M- It follows from (fT3T]) . (fT32|) and ffT26|) that 



= D2(f [X] + A" (°)/' v[X]) . (133) 

For future use we choose X = Xv in ( 1133^ and use the fact that v[Xv] = — Q, as 
follows from (1127^ . Equation fll33p assumes the form 

= (f [X,] - A" (°)/' Q) . (134) 
Relation between the variables 

We need an expression for the gauge-invariant linear perturbation H[X] of the Hub- 
ble scalar H of the preferred congruence, which is defined by 

H[X] = a((i)i7 -(°)i7'X0), (135) 

in accordance with the general definition ([3]). It follows from the expression (B.41a) 
for a^^'^H in UW, in conjunction with and f lT28D . that0 



3 

We have also used the fact that 



H[X] = iD2(v[X] - B[X]) - {d,^X] + n^X]). (136) 



i^''^ = 'H' -n^ = ~i^AG-K), (137) 
the second equality following from fl36l) . 



"^We note that "VaA^e) = °Va^(e) for a scalar A, and that °VjA = and "VqA = 
'*'*This is another example of a Replacement Principle, where an equation remains valid when each 
gauge- variant quantity is replaced by the associated gauge invariant defined using X-compensation. 
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First we note that 

v[X]-B[X] =V[X]-B[X]-Q = V-Q, (138) 

as follows from fll27p and the X-independent gauge invariant [V, B] in (1291) . Second 
we can use the transition rules (128|) and (132|) to show that the gauge invariant 

9^^[X] + H^X] + {-W + •H^)V[X], (139) 

is X-independent. Evaluating this expression for X° = X° and X° = Xp yields 

dr,^, + V,^, = -KY, (140) 

where we have used the linearized Einstein equation ( I40dp . the background Einstein 
equation At = Ac and the definition (136|) of Ac- Finally we choose X° = X° 
in ( I136P and substitute from fll38p and fll40p to obtain the desired result that 

3H^ = (D2 + 3K)V-D2Q = Z-D2Q, (141) 

the second equality following from ( HSj) . 

We can now use f ll34p to relate the perturbations of the variables D and Z in the 
1 + 3 approach to the corresponding variables D and Z in the metric-based approach. 
First choose f = p, X = M and n = 2 in ffTMp . and use ffTMD . ffTTTD . ffTOal) and fl35bll 
to obtain 

= D2(ro + 3'HQ). (142) 

Second, choose f = H , X = a and n = 1 in (11340 . Equations fll4ip and (I137p then 
lead tcS 

= (Z - (D^ + 3K~ ^Ag)Q) ■ (143) 

We also need the relation 

(^)g = D^Q, (144) 

which can be derived from the definition (1881) of Q and the definition f ll26p of Q. 
The desired equations f l98ap and ( ]98bp now follow immediately from f ll42p and f ll43p 
in conjunction with (I144p . Further, equations (I98cp can be derived from the defini- 
tions (190|) and (IHHj) of F and 11. Finally, equation (!97|) can be derived in a similar 
manner using (!96|) and the footnote following (I130p . 
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